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PREFACE 


This  paper  was  prepared  for  presentation  during  Session  No.  36  on 
Mechanics  of  Granular  Media  and  the  Related  Contact  Problems  at  the 
American  Society  of  Civil  Engineers  Annual  Convention  and  Exposition  in 
Atlanta,  Georgia,  22-26  October  1979. 

The  paper  describes  work  performed  by  personnel  of  the  Geomechanics 
Division,  Structures  Laboratory  (SL),  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  for  the  Defense  Nuclear  Agency  (DNA)  under 
Subtask  Y99QAXSB209,  "Propagation  of  Ground  Shock  Through  Soils  and 
Rock."  It  was  prepared  by  Dr.  J.  G.  Jackson,  Jr.,  Mr.  J.  Q.  Ehrgott, 
and  Dr.  Behzad  Rohani;  Dr.  Jackson  made  the  oral  presentation. 

Mr.  Bryant  Mather  was  Chief  of  SL  during  the  preparation  of  this 
paper.  The  Commander  and  Director  of  WE S  was  COL  Nelson  P.  Conover,  CE, 
and  the  Technical  Director  was  Mr.  F.  R.  Brown.  Dr.  George  W.  Ullrich 
was  the  Subtask  Manager  for  DNA. 
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INTRODUCTION 


Soil  conroressibility  is  usually  determined  in  the  laboratory  by  axially 

ii 

compressing  cylindrical  soil  samples  in  a  uniaxial  strain  device  and  measuring 
their  axial  deformation.  For  practical  applications,  the  laboratory  tests 
should  reflect  the  stress  levels,  loading  rates,  and  drainage  conditions  of 
interest. 

As  part  of  a  pioneering  research  effort  at  MIT  to  define  the  response  of 
soils  for  problems  associated  with  ground  shock  from  nuclear  explosions. 

Whitman  developed  a  dynamic  uniaxial  strain  device  based  on  a  novel 


1  Supv.  Res  Civil  Engr,  Geomechanics  Div,  Structures  Lab,  U.  S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  Mississippi. 

^  Res  Civil  Engr,  Geomechanics  Div,  Structures  Lab,  U.  S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  Mississippi. 

3  Res  Civil  Engr,  Geomechanics  Div,  Structures  Lab,  U.  S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  Mississippi. 

**  Although  frequently  referred  to  in  conventional  soil  mechanics  practice  as 
a  consolidometer ,  oedometer,  or  one-dimensional  compression  device,  the  test 
boundary  conditions  are  most  accurately  described  from  a  continuum  mechanics 
viewpoint  by  the  term  uniaxial  strain. 
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multiple-reflection  technique  (l,  2).  If  the  rise  time  of  the  applied  loading 
pulse  is  long  relative  to  the  time  required  for  a  stress  wave  to  propagate  back, 
and  forth  between  the  rigid-bottom  boundary  and  the  free-surface  boundary,  then 
inertial  stresses  can  be  neglected  and  stress  and  strain  within  the  specimen 
can  be  readily  deduced  as  a  function  of  time  from  external  measurements  of 
applied  pressure  and  surface  deflection.  The  MIT  gas-fluid  loading  system 
permitted  dynamic  pressures  up  to  about  1.1+  MPa  to  be  applied  with  rise  times 
on  the  order  of  15  msec.  Schindler  subsequently  developed  a  similar  test 
device  at  WES  which  used  a  piston- fluid  loading  system  to  develop  2.1-MPa 
pressures  in  about  3  msec  (3). 

Experimental  results  for  dry  sands  show  that  the  stress-strain  curves, 
starting  from  some  initial  prestress,  are  S-shaped,  with  yielding  for  small 
stress  changes  and  stiffening  for  large  stress  changes  (1»).  The  secant  modulus 
of  dry  sand  for  a  rapid  loading  (say  10-msec  rise  time)  is  usually  on  the  order 
of  5  to  1+0  percent  greater  than  that  for  a  slow  loading  (several  minutes  rise 
time);  typical  values  previously  reported  for  such  dynamic-to-static  modulus 
ratios  are  given  in  Table  1.  In  summarizing  the  MIT  findings.  Whitman  (8) 
concludes  that  time-dependent  effects  can  generally  be  ignored  for  problems 
involving  dry  granular  materials  subjected  to  dynamic  loadings  with  millisecond- 
range  rise  times,  provided  that  the  stress-strain  relation  used  in  the  analysis 
is  obtained  from  a  test  in  which  loading  and  unloading  take  place  in  30  msec  or 
less.  But  he  warns  that  there  is  evidence  to  indicate  that  time  effects  become 
very  important  when  the  duration  of  the  stress  pulse  drops  to  about  1  msec  or 
less  and  cites  a  need  for  more  adequate  study  of  time- dependent  effects  during 
uniaxial  strain  tests  with  large  stress  changes  and  submillisecond  rise  times. 


TABLE  1 — Secant  Constrained  Modulus  Ratios  Reported  for  Dry  Sands 


Sand  type  and 
source  of  data 
(1) 

Relative 

density 

% 

(2) 

Initial 

prestress 

MPa 

(3) 

Live 

stress 

MPa 

(*+) 

Dyn  stress 
rise  time 

msec 

(5) 

Byn-to-static 

modulus 

ratio 

(6) 

97 

0.10 

0. 3*» 

-10 

1.0*+ 

Uniform  beach  sand,  SP, 

38 

0.10 

0.3*+ 

-10 

1.10 

median  diam  =  0.3*+  mm 

70 

0.10 

0.3*+ 

-10 

1.15 

(Moore — ref  5) 

70 

0.10 

0.69 

8 

1.08 

70 

0.69 

0.07 

10 

1.38 

Uniform  fine  sand,  SP, 

median  diam  =  0.2*+  mm 
(Schindler — ref  6) 

9h 

0.0 

2.07 

25 

1.08 

Fine  silty  sand,  SM, 

median  diam  =  0.1*+  mm 

58 

0.0 

2.07 

5 

1.20 

(Jackson — ref  7) 


SCOPE  OF  THIS  PAPER 

Equipment  is  described  for  conducting  multiple-reflection  type  uniaxial 
strain  tests  to  much  higher  pressures  and  at  much  faster  loading  rates  than 
those  discussed  above.  Experimental  results  are  presented  from  tests  conducted 
on  three  dry  remolded  sands  (see  Fig.  l)  which  show  that  the  resistance  to 
compression,  or  stiffness,  of  these  materials  increases  dramatically  when  the 
loading  times  drop  below  a  millisecond.  Results  for  20-U0  Ottawa  Sand  are 
analyzed  with  a  one-dimensional  plane  wave  propagation  computer  code  to  deter¬ 
mine  if  the  measurements  were  adversely  affected  by  inertia  effects.  Results 
from  Uo  tests  on  HH  Clayey  Sand  are  analyzed  to  determine  if  the  observed 
trends  are  statistically  significant.  Finally,  results  from  tests  on  FH2 
Backfill  Sand  spanning  a  wide  range  of  submillisecond  loadings  are  compared  with 


the  response  of  a  three-element  rheologic  model. 


TEST  DEVICES  AND  TECHNIQUES 

The  100-MPa,  ram-loaded  uniaxial  strain  test  device  shown  in  Fig.  2  was 
developed  in  1969  and  is  similar  in  concept  to  the  device  described  in  detail 
by  Schindler  (6).  The  gas-driven  ram  can  apply  piston  loads  up  to  150,000  N 
with  rise  times  as  fast  as  3  msec  and  decay  times  as  fast  as  20  msec;  the 
piston,  in  turn,  produces  a  planar  pressure  pulse  in  an  oil-filled  chamber 
which  uniformly  loads  the  surface  of  the  membrane-sealed  soil  specimen.  Oil 
pressures  are  measured  with  a  flush-mounted,  diaphragm- type  pressure  trans¬ 
ducer.  The  specimens  are  laterally  constrained  by  the  thick  steel  walls  of  the 
bottom  assembly.  The  specimens  are  quite  thin  (1.27  cm  high),  in  order  to 
minimize  the  transit  time  of  propagating  stress  waves,  and  have  a  small  height- 
to-diameter  ratio  (1:7.6),  in  order  to  prevent  sidewall  friction  from  influencing 
center  deflections.  Center  deflections  are  measured  with  an  LVDT  transducer 
located  totally  within  the  pressurized  oil  chamber  and  mounted  directly  on  the 
soil  specimen  container. 


Fig.  2 — 100  MPa  Ram-loaded  Uniaxial  Strain  Test  Device. 


The  100-MPa,  explosive-loaded  uniaxial  strain  test  device  shown  in 
Fig.  3  was  developed  in  1978*  Instead  of  a  gas-driven  ram,  an  explosive 
charge  is  used  to  drive  a  piston  into  the  oil,  which  in  turn  pressurizes  the 
specimen.  The  device  consists  of  three  separate  assemblies  which  are  joined 
by  16  high-strength  bolts.  The  top  assembly  contains  the  explosive  charge 
and  chamber,  three  gas  expansion  chambers,  flow  controls  and  metering  orifices 
and  the  oil-pressurizing  piston.  Pressures  can  be  varied  up  to  100  MPa  by 
varying  the  size  of  the  charge;  rise  and  decay  times  as  fast  as  0.3  msec  can 
be  obtained  with  the  gas  flow  controls.  The  middle  assembly  contains  the  oil 
chamber,  the  flush-mounted  pressure  transducer,  and  the  electrical  leads  to 
the  LVDT  transducer.  The  lower  assembly  holds  the  test  specimen  container. 

The  same  soil  specimen  container  and  deflection  measurement  system  were 
used  in  both  the  explosive-loaded  and  the  ram-loaded  devices.  Preweighed 
quantities  of  sand  were  alternately  spooned  and  tamped  into  the  container  and 
the  surface  leveled  with  a  straightedge.  A  5-mm-thick  rubber  membrane 
attached  to  an  aluminum  LVDT  footing  assembly  was  placed  over  the  specimen 
and  sealed  around  the  edges  with  a  steel  clamp  containing  two  0-rings. 
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TESTS  OH  20-10  OTTAWA  SAND 

Six  tests  were  conducted  on  air-dried  specimens  of  20-t0  Ottawa  Sand  with 
an  initial  void  ratio  of  0.572;  stress-strain  and  pressure-time  plots  are 
given  in  Fig.  5.  Results  from  the  two  static  tests  are  in  excellent  agreement 
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with  each  other  (see  Fig.  5a);  they  also  agree  quite  well  with  previously 
reported  high-pressure  data  for  a  20-1*0  Ottawa  Sand  with  an  initial  void  ratio 
of  0.527  (9) •  Results  from  two  dynamic  tests  with  the  ram-loaded  device  are 
shown  in  Fig.  5b;  the  specimen  in  test  D1  was  loaded  to  30  MPa  in  about 
100  msec  while  that  in  test  D2  was  loaded  to  Uo  MPa  in  about  1*  msec.  The 
stress-strain  curves  for  D1  and  D2  are  essentially  identical;  they  could  also 
be  easily  mistaken  for  the  SI  and  S2  results.  The  average  dynamic-to-static 
ratio  of  secant  modulus  to  30  MPa  computed  from  these  data  is  1.06,  a  value 
quite  consistent  with  those  given  in  Table  1  for  low-pressure  tests.  But  the 
results  shown  in  Fig.  5c  for  dynamic  tests  D3  and  DU,  in  which  1+0-MPa  pressures 
were  applied  in  0.3  to  0.1*  msec  with  the  explosive-loaded  device,  are  markedly 
different,  both  qualitatively  and  quantitatively,  from  the  other  results. 

Since  the  frequency  response  of  the  measurement  system  was  above  15  kHz,  the 
measurements  should  be  valid.  However,  before  this  difference  in  response  can 
be  attributed  to  loading  rate  effects,  an  analysis  is  necessary  to  insure  that 
it  was  not  due  to  wave  propagation  or  inertia  effects. 

INERTIA  ANALYSIS 

A  computer  code  for  analyzing  one-dimensional  plane  wave  propagation 
through  layered  profiles  of  nonlinear  hysteretic  materials  (10)  is  also  used  at 
WES  to  assess  inertia  effects  in  multiple-reflection  uniaxial  strain  tests. 

The  method  of  analysis  is  based  upon  a  discrete  model  consisting  of  lumped 
masses  connected  by  springs  and  subjected  to  an  arbitrary  dynamic  pressure  at 
the  free  surface  (Fig.  6).  Piece-wise  linear  stress-strain  relations  for  both 
loading  and  unloading  are  used  to  define  the  spring  resistances. 
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MEASURED  RESULTS 


FOR  OTTAWA  SAND  TESTS 
H  =  1.27  cm 
A  =  73.1  cm2 
7  =  1.69  gm/cm3 

RESULTS  FROM 
STATIC  TEST  S2 

EISED  TO  define] 

1  FOR  ID  CODE 

:alculation  J 


CALCULATED  o-e  RELATION 
USING  0.5  iiSEC  TIME  STEP 
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o  =  Pit) 
uJt) 

e=-9TX,0° 


STRAIN,  % 


Fig.  6 — Lumped-parameter  Fig.  7 — Inertia  Analysis  for  Dynamic  Test  DU  on 

Model  for  Analyzing  20-^0  Ottawa  Sand  Assuming  Rate-independent 

Multiple-reflection  (Static)  Stress-strain  Properties. 

Uniaxial  Strain  Tests. 

Inertia  influences  within  the  20-1*0  Ottawa  Sand  tests  can  be  bounded  by 
using  the  static  test  results  to  define  rate-independent  properties  for  the 
springs,  and  then  subjecting  the  system  to  the  most  severe  dynamic  pressure 
pulse.  The  analysis  performed  for  dynamic  test  Dl*  is  shown  in  Fig.  7.  The 
differences  between  the  o-e  relation  based  on  surface  pressure  and  cal¬ 
culated  near-surface  displacements  and  the  o-e  relation  from  static  test 
S2  (which  was  specified  as  the  "true"  property)  represent  the  maximum  potential 
measurement  errors  due  to  inertia  effects.  Thus  the  measured  results  from 
tests  D3  and  D4  on  20-h0  Ottawa  Sand  were  not  significantly  affected  by  wave 
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propagation  within  the  specimen;  they  apparently  depict  a  real,  and  rather 
dramatic,  change  in  material  behavior  as  a  result  of  the  submillisecond  loading 
times. 

TESTS  ON  HH  CLAYEY  SAND 

Forty  tests  were  conducted  on  densely  compacted  specimens  of  HH  Clayey 
Sand;  nineteen  of  the  specimens  were  statically  loaded  to  10  MPa  in  1-10  min, 
sixteen  were  ram-loaded  to  10  MPa  in  5-100  msec,  and  five  were  explosively 
loaded.  Peak  pressures  applied  during  the  latter  tests  ranged  from  ko  to 
6o  MPa  with  the  initial  load-unload  cycle  completed  in  slightly  less  than 
1  msec;  results  are  shown  in  the  left  half  of  Fig.  8. 


Fig.  8— Uniaxial  Strain  Test  Results  for  HH  Clayey  Sand. 
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Although  the  results  appeared  to  indicate  appreciable  rate  effects,  they 
were  obscured  somewhat  by  the  data  scatter  inevitably  present  in  compacted 
specimens  containing  significant  fractions  of  silt  and  clay-size  particles.  In 
order  to  isolate  the  rate-effect  trends ,  the  strain  recorded  in  each  test  at  a 
stress  of  1  MPa  was  subtracted  from  the  total  strain  recorded  at  various  stress 
levels  up  to  10  MPa.  Mean  values  and  standard  deviations  were  computed  for 
each  of  the  three  test  groups;  these  are  plotted  in  the  right  half  of  Fig.  8. 
Analysis  indicates  that  the  groupings  reflect  three  statistically-distinct 
populations.  Dynamic-to-static  secant  moduli  ratios  computed  from  the  mean 
strain  values  at  10  MPa  are  1.1*5  for  the  5-100  msec  rise  times  and  a  whopping 
6.28  for  the  0.2-0. 3  msec  loadings.  Tests  need  to  be  conducted  at  several 
intermediate  loading  rates  in  order  to  quantify  the  transition  between  these 
values . 

TESTS  ON  FH2  BACKFILL  SAMP 

Fifteen  tests  were  conducted  on  FH2  Backfill  Sand  (see  Fig.  9)*  The 
slowest  loadings  were  applied  in  test  SI  (about  20  MPa  in  35  sec)  and  test  S2 
(about  1+0  MPa  in  50  sec).  Tests  S3  and  Sl+  were  performed  with  loadings  to  35 
MPa  in  about  20  sec.  Results  for  these  four  static  tests  are  plotted  in  Fig. 
9a.  The  specimen  in  dynamic  test  D1  was  subjected  to  an  80-MPa  pressure  in 
about  1+  msec  with  the  ram  loader.  By  varying  the  size  of  the  charge,  the 
weight  and  size  of  the  piston,  and  the  flow  control  settings,  it  was  possible 
to  produce  a  range  of  loading  rates  with  the  explosive- loaded  device,  e.g. ,  10 
MPa  in  about  1  msec  in  tests  D2  and  D3,  15  MPa  in  about  0.5  msec  in  tests  Dh 
and  D5,  and  30  MPa  in  0.5  msec  in  tests  d6  and  D7.  Results  for  these  seven 
dynamic  tests  are  plotted  in  Fig.  9b.  Finally,  dynamic  tests  D8-D11  were 
performed  with  pressure  pulses  which  reached  30  MPa  in  about  0.3  msec;  these 
results  are  plotted  in  Fig.  9c. 


The  strain  recorded  at  10  MPa  was  again  chosen  as  a  basis  for  quantifying 


loading  rate  influences  on  compressibility.  The  static-to-dynamic  strain 
ratio  at  10  MPa  for  each  test  is  plotted  versus  loading  time  in  Fig.  10.  The 
plot  demonstrates  conclusively  that,  while  rate  effects  on  the  compressibility 
of  this  dry  sand  can  generally  be  ignored  for  loading  times  greater  than 
1  msec,  they  cannot  be  ignored  for  submillisecond  loadings,  i.e.,  stiffness 
increases  rapidly  as  loading  time  decreases  within  the  1-msec  to  0.1-msec 
decade.  Thus,  the  suspicions  of  Whitman  (8)  appear  to  have  been  well-founded. 

The  dynamic-to-static  modulus  ratios  at  30  MPa  are  not  nearly  so  extreme 
(a  maximum  of  about  1.5) »  but  the  overall  stress-strain  response  shows  a 
marked  change  in  character  as  the  loading  rate  is  increased  above  about 


Fig.  10 — Effect  of  Loading  Rate  on  Compressibility  of  FH2  Backfill  Sand. 


10  MPa/msec.  Belov  this  rate,  the  stress-strain  curve  "stiffens"  (curvature 


whereas  above  this  rate,  the  curve 


Such  changes  in  curvature  would  he're 


significant  effects  on  propagation  velocities  and  rise  times  at  the  fronts  of 


blast-induced  stress  waves  ( k ,  11).  A  variety  of  simple  viscoelastic  and 


hysteretic  models  have  been  used  to  study  the  effects  of  time-dependent 


inelastic  and  nonlinear  stress-strain  behavior  on  stress  wave  propagation 


through  laterally  confined  soils  (12,  13);  a  model  incorporating  all  of  these 


characteristics  is  needed  in  order  to  replicate  the  behavior  shown  in  Fig.  9 


THREE-ELEMENT  RHEOLOGIC  MODEL  CALCULATIONS 


The  three-element  rheologic  model  used  in  this  study  is  a  modified  version 


of  a  three-element  viscoelastic  model  described  by  Kolsky  (lU)  in  which  a 


linear  spring,  M  ,  is  coupled  in  series  with  another  linear  spring,  M  ,  in 


parallel  with  a  dashpot,  n.  For  static  loadings,  the  response  of  such  a 


The  WES  modification  involved  a 


nonlinear-hysteretic  formulation  of  M. 


so  that  M 


the  complete  stress-strain  curve  measured  during  the  static  tests  after  a 


value  for  M0  has  been  fixed  to  represent  the  stiff  "elastic"  response 


depicted  by  unloading  data.  The  governing  differential  equation  of  the  model  is 


solved  numerically  for  an  arbitrary  stress-time  input  using  values  of  M. 


Whenever  the  total  strain  increment  is 


positive,  these  values  are  calculated  from  the  loading  portion  of  the  Mf 


relation;  otherwise  they  are  calculated  from  the  unloading  portion 


The  model  is  pictured  on  the  right-hand  side  of  Fig.  11a  along  with  the 


spring  and  dashpot  coefficient  values  used  in  calculating  its  response  to  the 


pressure  pulses  measured  during  the  FH2  Backfill  Sand  tests.  The  2000-MPa 


value  for  M  was  obtained  from  the  static  unloading  data;  M  was  derived 


.U 


-y.  "•!  > 


from  the  piece-wise  linear  hysteretic  fit  for  ,  which  is  shown  on  the  left 
for  comparison  with  results  from  static  tests  SI  and  S3.  The  value  of 
5/1^5  MPa-sec  for  n  was  selected  after  a  series  of  parametric  (or  trial-and- 
error)  calculations .  The  model  was  then  driven  by  three  loading  functions 
based  on  the  pressure  pulses  measured  during  dynamic  tests  D2-D7.  These  func¬ 
tions  are  depicted  by  the  heavy  solid  lines  in  Fig.  lib;  the  measured  pressure¬ 
time  histories  are  shown  for  comparison.  The  calculated  stress-strain 
responses  are  also  shown  as  heavy  solid  lines  in  Fig.  lib  for  comparison  with 
the  measured  results.  The  agreement,  both  quantitatively  and  qualitatively,  is 
very  good.  The  stress-strain  curves  from  tests  D8-D11  exhibit  a  very  erratic 
character  consisting  of  alternating  "bumps"  and  "flats"  (see  Fig.  9c)  which 
were  thought  at  first  to  represent  either  measurement  inaccuracies  or  wave 
propagation  effects.  However,  when  the  model  is  driven  with  the  pressure-time 
pulses  actually  measured  in  tests  D8  and  DIO  ( as  shown  in  Fig .  11c ) ,  the 
calculated  stress-strain  curves  reflect  the  detailed  character  of  the  measured 
results  remarkably  well.  It  thus  appears  that  the  observed  "bumps"  and  "flats" 
are  not  due  to  measurement  inaccuracies ,  but  rather  reflect  real  rate-dependent 
behavior. 

SUMMARY  AND  CONCLUSIONS 

Experimental  results  are  presented  on  the  uniaxial  strain  response  of 
three  dry  sand  materials  to  intense  transient  pressure  pulses;  loadings  ranging 
from  10  MPa  to  SO  MPa  were  applied  with  rise  times  ranging  from  a  few  tenths  of 
a  millisecond  to  a  few  minutes.  These  data  are  consistent  with  previous 
findings  regarding  the  dynamic  compressibility  of  dry  sand,  in  that  loading 
rate  had  a  relatively  minor  effect  on  the  stress-strain  response  for  loading 
times  greater  than  1  msec.  They  also  confirm  previous  suspicions  that  this 
would  not  be  the  case  for  submillisecond  loadings,  in  that  secant  modulus  values 


at  10  MPa  increased  by  an  order  of  magnitude  within  the  1-msec  to  0. 1-msec 
rise-time  decade. 

In  addition  to  the  dramatic  quantitative  effect,  the  overall  response 
depicted  a  marked  change  in  character  as  the  loading  rate  was  increased  above 
10  MPa/msec,  i.e.,  stress-strain  curvature  rapidly  shifted  from  "stiffening" 
(concave  to  the  stress  axis)  to  "yielding"  (concave  to  the  strain  axis). 
Stress-strain  curves  obtained  with  a  relatively  simple  rate-dependent  rheologic 
model  reflect  the  detailed  character  of  the  measured  results  remarkably  well. 
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